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Abstract 

We pursue the issue of the local quark-hadron duality at high energies in two- 
and four-dimensional QCD. A mechanism of the dynamical realization of the 
quark-hadron duality in two-dimensional QCD in the limit of large number of 
colors, N c — > oo, (the ' t Hooft model) is investigated. We argue that a similar 
mechanism of dynamical smearing may be relevant in four-dimensional QCD. 
Although particular details of our results are model-dependent (especially in 
the latter case), the general features of the duality implementation conjectured 
previously get further support. 

I. INTRODUCTION 

In the recent years the focus of applications of the operator product expansion (OPE) 
has shifted towards the processes with the essentailly Minkowskean kinematics. Perhaps, 
the most well-known example is the theory of the inclusive decays of heavy flavors (for a 
review see e.g. Ref. ||). This fact, as well as the increasingly higher requirements to the 
accuracy of predictions, puts forward the study of the quark— hadron duality as an urgent 
task. 

A detailed definition of the procedure which goes under the name of the quark-hadron 
duality (a key element of every calculation referring to Minkowskean quantities) was given 
in Refs. [0,|J. 111 a nu t shell, a truncated OPE is analytically continued, term by term, from 
the Euclidean to the Minkowski domain Q A smooth quark curve obtained in this way 
is postulated to coincide at high energies (energy releases) with the actual hadronic cross 
section. 

If duality is formulated in this way, it is perfectly obvious that deviations from duality 
must exist. In Ref. || it was shown that one of the sources is the asymptotic divergence of 
high orders in the power series. If we knew the leading asymptotic behavior of the high order 
terms in the power series we could predict the pattern of the duality-violating contributions 
at high energies. Unfortunately, very little is known about this aspect of OPE, and we have 
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to approach the problem from the other side - either by modelling the phenomenon |4[] or 
by studying some general features of the appropriate spectral densities. One can also try to 
approach the problem purely phenomenologically (for recent attempts see e.g. |5|-0). 

An illustrative spectral density, quite instructive in the studies of the issue of the quark- 
hadron duality, was suggested in ||, 

oo 

ImIT = Const.iV c tt ^ S(£ - n) (1.1) 

n=l 

where 

and we droped an inessential constant in front of the sum. The color factor N c is singled 
out for convenience. The imaginary part above represents, for positive values of s, a sum of 
infinitely narrow equidistant resonances, with equal residues. It defines II(g 2 ) everywhere in 
the complex plane q 2 , through the standard dispersion relation, up to an additive constant 
which can be adjusted arbitrarily. It is not difficult to see that the corresponding correlation 
function 



II(g 2 ) = -N c 



;i-2) 



where ip is the logarithmic derivative of Euler's gamma function, and 

In the Minkowski domain £ is positive, in the Euclidean domain e is positive. Then, the 
asymptotic expansion of II(g 2 ) in deep Euclidean domain is well known, 



IT(g 2 ) - —N c 



1 B 2n e 



-2n 



2e 2n 



(1.3) 



where B 2n are the Bernoulli numbers. At large n they grow factorially, as B 2n ~ (2n)\ (see 
||, page 23), and are sign alternating. 

Although the spectral density ( |1 . 1|) is admittedly a model, it was argued || that a similar 
factorial growth of the coefficients in the power (condensate) series is a general feature. 

The spectral density ( |1 . 1| ) may be relevant in the limit of the large number of colors, 
N c — > oo, when all mesons are infinitely narrow. This limit is not realistic, however, and, 
moreover, in this limit the local quark-hadron duality, as we defined it, never takes place 
since even at high energies the hadronic spectral density never becomes smooth. One can 
smear it by hand, of course, but then deviations from the local duality will be determined 
not only by the intrinsic hadronic dynamics, as is the case in the real world, but also by 
details of the smearing procedure - the weight function chosen for smearing, the interval 
of smearing and so on. In the actual world the smearing occurs dynamically, since at high 
energies the resonance widths become non-negligible. The limits of E — > oo and N c — > oo 
are not interchangeable. 
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In this work we study more realistic (dynamically smeared) spectral densities compatible 
with all general properties of Quantum Chromodynamics. Starting from infinitely narrow 
resonances, as in Eq. fll.lj) , we then introduce finite widths, ensuring smooth behavior. 
Various dynamical regimes leading to specific patterns of duality violations are considered, 
with a special emphasis on the high-order asymptotics of the power series. Technically, 
in the first part of the paper the problem of duality is analyzed in the two-dimensional 't 
Hooft model || (see also ]TT|-|T3"||). We see that there indeed exists the dynamic smearing 
of the spectral density. This smearing occurs already if one takes into account the fact 
that the resonance widths are nonzero and calculates the relevant spectral densities in the 
Breit-Wigner approximation. We see that there are two characteristic scales in the problem. 
Starting from the first scale the spectral density is well approximated by the oscillations 
round the average value given by the OPE asymptotics. Starting from the second scale the 
oscillations become small and local duality is valid. 

The ^-function model was originally suggested in Ref. || for the heavy-light quark 
systems. In Ref. |14| it was noted that it was more appropriate for the light-light quark 



systems since in the heavy-light systems the resonances are not expected to be equidistant. 
Straitforward quasiclassical estimates yield in this case that the meson energies (measured 
from the heavy quark mass) asymptotically scale as yfn. (More elaborated models, e.g. the 
one of Ref. | 15[| , based on a relativistic linear potential, yield not only the y/n law, but also 



the coefficient in front of yjn.) In the present paper we further develop the ?/>-function model 
adapting it for the light-quark systems. 



II. RESONANCE WIDTHS IN T'HOOFT MODEL 

Our main original contribution to t 'Hooft model is the calculation of the first nonvanish- 
ing corrections to the widths of the resonances of t'Hooft model. Then, using these widths, 
we shall find the asymptotic behavior of the polarization operator of two scalar currents. 
This resonance-saturated polarization operator will be referred to as phenomenological. We 
will confront it with the truncated power expansion. The difference between these two 
expressions presents duality violation. 

It was shown by t'Hooft that the 2d QCD is exactly solvable in the limit N c — > oo. 
The bound state spectrum includes an infinite number of bound states whose masses m n lie 
on the almost linear trajectory. The properties of these bound states are described by the 
t'Hooft equation 

2 / / \ (7 2 - i) f 1 Mv)dy fou 

Here the integral is understood as a principal value, The variable x is the momentum fraction 
of the meson carried by the antiquark (in the infinite momentum frame), while 1 — x is that 
of the quark; 4> n (x) is the wave function of the n-th bound state, [i 2 n = m^/fi 2 



The integral equation |2.1| must be solved with the boundary conditions: 

4> n {x) -> X P(l -xf (2.2) 

for x — > 0, 1. Here 
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7r/3ctg(7r/3) = 1 - 7 2 (2.3) 

and 7 2 = m 2 //i 2 , m q is the mass of the quark. 

Below we shall be interested in the massless case 7 = and in the number of flavors 
equal to 1. In this case the asymptotic behavior of /z 2 versus n is given by equation 

/i 2 =ir 2 n(l+0(\n(n)/n + ...) . (2.4) 

As for the wave functions <f> n (x) , their calculation is quite complicated even numerically, 
especially in the massless case. The recent calculations ]T7| , exploiting a new improved 
numerical procedure (spline method), show that for the massless case these wave functions 
are very accurately approximated by 

4>n(x) = \Z2cos(irnx) . (2.5) 



Note that the wave functions |275| satisfy the proper boundary conditions for the massless 
case. 

Let us now calculate the meson widths. As it was already mentioned, in the limit 
N c — > 00 the bound states in the t'Hooft model are stable, their widths vanish. However, 
once one takes into account the leading 1/N C correction, the resonances begin to decay. In 
the first order in 1/N C expansion there are only two-particle decays a —>■ b + c. The relevant 
coupling constatnts g abc are given by jll],[n^[n].[l^,[n| : 



9abc 



{ -^^/N c (l - (-l)^ + ^)(fL + f~ bc ) (2.6) 

7T V 



Here a a is the parity of the a-th resonance. The constants f^ bc are determined from the 
following expressions: 

1 r UJ ± I x jj.. \ 

fabc = ~, / (j) a {x)(j) b {x/u} ± )^ c dx , 

I — LU± JO \ 1 — Ul± I 



— ^ a (x)$ b (x/w ± )0 c f dx . . (2.7) 



Here u± are two roots of the algebraic equation corresponding to the mass-shell condition, 

ml = ^ + -^- y (2.8) 
uj {1 — UJ) 

Two different idices of uj correspond to two different modes of the two-dimensional decay: 
in the rest frame of the resonance a the resonance b can go to the right and the resonance c 
to the left, and vice versa. 

The function $ a (^) is defined as 

$ a (x) = f dy(j) a {y)/{x - yf. 
Jo 
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It is easy to carry the calculation of the couplings f a bc analytically (see Ref. jlj) and ex- 
press them via integral sinus and integral cosinus functions. Using the above expresions 
for the decay couplings one can readily calculate the resonance widths in the leading 1/N C 
approximation. They are given by 



(m 2 a + m 2 b - m 2 c )g, 



Y _ _ \""a 1 ""ft ""cJUabc ^ 9) 



abc \ i(mi,mt,mi 



where I is the standard "triangular" function, 

1 



I(x,y,z) = ^[m 2 a - (m b + m c ) 2 ][m 2 a - (m b - m c ) 2 } . (2.10) 



The sum in Eq. fl2.9|) runs over all mesons b and c whose masses are lighter than that of a. 

Details of our calculation of the resonance widths are described in Ref. 0. The calcula- 
tion was done numerically using the trial wave functions (|2.5| ). Our task was to establish the 



asymptotic behavior of the widths, as a function of the excitation number, at large values of 
n, in the leading 1/N C approximation (all widths are proportional to 1/N C ). After computing 
the overlap integral we performed the summation over b and c. The result for the widths 
exhibits a remarkable pattern. The widths of the individual levels oscillate near a smooth 
square- root curve, see Fig. 1. This plot shows the width of the a-th state versus a, up to 
a = 500. The result of averaging over the interval of 20 resonances is depicted in Fig. 2. 
We see that the curve of the resonance widths T(a) = r a is very well approximated by the 
function 

' 1/4 y/E(l + 0(l/a)) , (2.11) 



ir 3 N, 



where [i 2 = g 2 N c /n. 

Since the square-root law ( 2.11|) for the (averaged) widths is valid in such a large interval 



of the excitation numbers and turns out to be so accurate, it seems to be doubtless that this 
formula could be obtained analytically. This is an interesting question by itself, especially in 
the four-dimensional QCD. Unfortunately, we were unable to find analytic solution so far. 
The numerical value of the constant A is 

A ~ 0.44. (2.12) 

Below the above result for the (averaged) widths will be used for determining of the asymp- 
totic behavior of the polarization operator. 

Concluding this section let us note that the same square-root was reported previously in 
Ref. inj . We failed to reproduce the arguments of this work leading to the square- root law, 



however. Moreover, what is even more important, the constant A in Ref. |19| is claimed to 
be proportional to l/\/m (!), and, thus, blows up for massless quarks. This poses perplexing 
questions. The coincidence looks completely accidental. 



5 



III. QUARK-HADRON DUALITY IN 2D QCD 



Once we had found the resonance widths, we can calculate the polarisation operator in 
the Breit-Wigner approximation. Here we shall consider the most interesting case of the 
polarisation operator of the two scalar currents. Let us start from this polarisation operator 
in the N r — > oo limit HTT1,|T8 . 



The polarisation operator is given by the correlator 



dV* <0\T{j(x),j(0)}\0> 



(3.1) 



Here j(x) is the scalar current: 

j(x) = q(x)q(x) 

and q(x) is the quark field. It was shown in Ref. that this polarisation operator can be 



(3.2) 



calculated explicitly in the N c — > oo limit and is given by 

J n 



H(g* 



n=l 



mi + ie 



It was shown in Ref. JT8[ that in order for eqs. [372] and the perturbation theory for QCD 2 to 
be compatible, the coeffficients f n for sufficiently large n must be independent of n. Then, 
taking into account the linear dependence of the mass squared on n one can approximate 



the polarisation operator |3J] for sufficiently large q 2 by ^-function [[18 

H(g 2 ) - H(0) 



Nr. 



■It 2 / 2 2\ 

~ — ip{—q /7T /i 

71 



(3.3) 



f 2 ~ f 2 



H 2 TT 2 . 



Here lim r 

Consider now what will happen if we take into account the finite widths of the resonances. 
We shall calculate the polarisation operator in the Breit-Wigner approximation. In order to 
do it we shall first find the inverse propagator of the n-th bound state: 



n, 



+ ie) = q 2 -m 2 n + S(g 2 + ie) 



(3.4) 



Using the procedure described in Ref. [0 it is easy to prove (see Ref. |L] for details) that 



Z(q 2 + ie) 



Bq\ , q 2 
-Log(- — 



Nn 



ie 



(3.5) 



The latter equation is valid for q 2 > q 2 , where q 2 is some intermediate scale, and B = A/tt a . 
Then one obtains for the polarisation operator 



n(g 2 ) 



N 



71 



E 



/, 



f~< + £Log(- 



le) 



(3.6) 



(for the nonphysical sheet). Now we can use the known aymptotic behaviour of m 2 n and /, 
to sum the series [376] . We immediately obtain: 
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N c q 2 + ^Log(g 2 /g 2 ) + iBq 2 /N c 
I% 2 ) - n(0) ~ —^(-- 3* — -) (3.7) 

In other words we obtain that the polarisation operator on the unphysical sheet is the psi- 
function of the complex argument. This is our main result. In the limit N c — > oo we clearly 
recover eq. [O] . 

We can now study the duality violation and the operator product expansion for the 
polarisation operator |3.1| . 

Let us start from the imaginary part. The are two possible patterns of behaviour of 
the imaginary part of the polarisation operator as a function of q 2 : Euclidean kynematics 
corresponds to q 2 < and Minkowskean kinematics corresponds to q 2 > 0. Note that in 
the Eclidean domain the imaginary part of the polarisation operator is zero. Indeed, 



the logarithm in the argument of the ip function does not have any imaginary part, and the 
entire polarisation operator is real. Thus, the imaginary part of the polarisation operator is 
zero for q 2 < as it must be. Consider now the imaginary part in the Minkowskean domain 
(i.e. the spectral density). In the limit N c — > oo the imaginary part of the polarisation 
operator is the sum of the delta-functions, corresponding to the poles of the psi-function: 

f 2 N 

im(IT(g 2 ) - 11(0)) ~ — - V 5{q 2 - n/Ar 2 ). (3.8) 

TT „ 

However, once we take into account the leading 1/N C correction, the situation immediately 
changes. Instead of the sum of the delta functions, one obtains the oscillations of the 
imaginary part round the constant value N c /tt with decreasing amplitude. Note that the 
latter is just the average value of the imaginary part of |3]8| over the interval q 2 . Thus we 
see the example of dynamic smearing: instead of infinite peaks, QCD dynamically produces 
decreasing oscillations round average smooth function. 

Let us now check that this is indeed the case for the concrete example of our polarisation 
operator |3.6| . In order to find the relevant imaginary part explicitly, we shall make use of 
the reflection property of ^-function: 

1p(z) = Ip(—Z) — 7TCtg(7T2;) — 1/z. (3.9) 

Then, using eq. one obtains 



Im(n(g 2 ) - 11(0)) = ^{(Irm/H(g 2 + ^f- \og{q 2 /q 2 ) + ^/(ttV )) 
+ Im7rctg(vr(g 2 + ^- \ogq 2 /q 2 + iBq 2 / N c ) / (tt 2 fi 2 )) 

iV c 7T 

+ M „ — ^ )} 

q 2 + ^\ogq 2 /q 2 + iBq 2 /N c ' S 

(3.10) 

It is easy to see that a sum of the first and the third terms in eq. |3.10| is the smooth 
function of q 2 (for positive q 2 ). Numerically these functions are well approximated by the 
OPE (asymptotics). 
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Consider now the second term in eq. |3.10| . For N c 
the points where 



oo this term has simple poles in 



q 



1 1 

[inn, 



i.e. the imaginary part of this term is just the sum of the delta functions TB . 

Consider now the case of the large but finite N c . Then for q 2 > n 2 [i 2 the imaginary part 



of the polarisation operator |3.10| is well approximated by the imaginary part of the ctg term 
in eq. |3.10| : 



Im(n(g 2 ) - 11(0)) ~ -N c 



sh(2y) 



ch(2y) — cos(2a;) 



This means that for y — > oo 

i m n(g 2 ; 



-N c (l + 2exp (-2y)cos(2x) + ... 



(3.11) 



(3.12) 



Here x = q 2 /(n[i 2 ) + 0(1/ N c ) and y = Bq 2 /(N c ttli 2 ). 

We can now see the behaviour of the imaginary part of the polarisation operator. For 
q 2 > ix 2 [i 2 this behaviour is well approximated by the ctg term (eq. |3.11| ). The poles in the 
imaginary part of eq. [3.11 disappear for arbitrary nonzero y, i.e. for every finite N c . Instead 



we have peaks with the amplitude given for q ~ [i n by N c /(Bq ) The amplitude of these 
peaks goes to oo as N c — > oo. 
For 



q 2 > N c ttli 2 /(2B) 



(3.13) 



the asymptotic behaviour will be given by eq. [3.12| . The oscillations decrease exponentially 
and local duality is established. At these q 2 one may use the OPE. For q 2 smaller than 
the scale given by eq. |3.13j the oscillations dominate and we can speak only about global 
duality. For q 2 much larger than this scale one can safely use the OPE. The corresponding 
inaccuracy decreases exponentially as exp(—2q 2 B/(N c 7i[i 2 )). 

Numerically the typical behaviour of the spectral density is depicted in Fig. 2 for [i 2 = 
0.01, N c = 10. The function Fl depicts the (normalised) behaviour of the full ip function 
expression for the spectral density, F2-of the ctg term in eq. [3.7| and F3-the exponentially 
decreasing term F3. 

For q 2 < the polarisation operator is 



Consider now the full polarisation operator 3.6 



evidently given by the ip function of the positive real argument and is a smooth function 
with no singularities. It is well approximated by its asymptotic expansion 3.12: 



n(g 2 )-n(0) = ^(log(-g 2 /^ 

7T 



2n q 2 



(3.14) 



Here Bm are the Bernulli numbers. Note that the expansion [3.14 is not Borel summable. 



Indeed, for large n the ratio a n+ x/a n , where a n is the n-th term of the series [3.14j , becomes 
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~ 1 staring from n ~ — q 2 / (/i 2 7r) This means that the asymptotic expansion has the intrinsic 
uncertainty 



An(g 2 ) ~exp(^log-2g 2 /A 2 ). (3.15) 



Note that it is senseless to continue analytically the exponentially decreasing term |3.15| into 
the Minkowski domain, although it has the multiplier that begins to oscillate upon such 
continuation. This teaches us that it is extremely dangerous to continue different parts of 
the polarisation operator analytically from Euclidean to Minkowski domain. 

The 1/N C terms in the Euclidean domain are only small corrections to the N c — ► oo 
expression and can be neglected numerically. 

We can now carry the staightforward analysis of the large q 2 behaviour of the asymptotic 
expansion of the polarisation operator. This is easily done in the same way as it was done 
for the imaginary part. 

We immediately see that in the Minkowski domain the polarisation operator is given 
by the sum of the asymptotic expansion |3.14| (continued analytically into the Minkowski 



domain) and the oscillating part that for q 2 > N c irfi 2 / (2B) is equal to 

n osc (g 2 ) ~ 2— sin (2q 2 /(nfi 2 )) exp (-2Bq 2 /(N c 7r^ 2 )). (3.16) 

7T 

We once again see an oscillating part with the exponentially decreasing amplitude of the 
oscillations that must be added to the smooth OPE function. 

We conclude that there are indeed two characteristic scales in the problem. For q 2 > 7r 2 /z 2 
the polarisation operator is dominated by oscillations, and this oscillations begin to decrease 
exponentially at the scale given by eq. |3.13| . Starting from the latter scale the amplitude of 
the oscillations become smaller than the leading term in the OPE and there is sense to speak 
about local duality, that establishes itself for q 2 much bigger than the scale of eq. |3. 13 . The 



corresponding inaccuracies decrease exponentially. On other hand the global duality holds 
for all values of q 2 . 

IV. OPE AND THE POLARISATION OPERATOR IN 4D QCD. 

In this section we shall consider the OPE in 4d QCD. Once again, we consider the 
polarisation operator of 2 vector currents and calculate it for large absolute values of q 2 for 
both Euclidean and Minkowski domains. 

A. Polarisation operator of vector currents in the Veneziano model 

Our goal in theis subsection will be to study the polarisation operator of two vector 
currents: 

IV(g 2 ) = J exptqx < 0\T{j„(x)j v {0)}\Q > d A x. (4.1) 
Here j^{x) is the vector current, 
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jn(x) = uj^dix). (4.2) 
Due to the conservation of vector current the polarisation operator can be represented as 

n F (? 2 ) = (9V-9 2 v)n(g 2 ). (4.3) 

The polarisation operator can be calculated in the Euclidean domain using the OPE: 

n(g 2 ) ~log-g 2 + a /g 4 + .... (4.4) 

Usually in order to calculate II(g 2 ) in the physical, Minkowski domain we just analytically 
continue eq. [4.4| into the Minkowski domain. In particular, the logarithm acquires the 
constant imaginary part. Then we get 

n(g 2 ) ~7r + 0(l/g 2 )... (4.5) 

This is not the whole story however. Indeed, we know from our experience with QCD2 
that in addition to the smooth constant there may be oscillations. 

The question of oscillations is especially interesting since in QCD the imaginary part 
of the polarisation operator of vector currents is directly related to the cross-section of the 
e + e~ annihilation into hadrons. 



B. Polarisation operator and resonance widths 

We have seen already the connection between the resonance masses and widths asymp- 
totics and the asymptotic of the polarisation operator in the Breit-Wigner approximation in 
2d QCD. The same connection exists in 4d QCD. The problem is, however, that we do not 
know how to calculate the relevant resonance masses and widths in a model independent 
way, like it was done in QCD 2 . We need to know the masses and widths of the radial exci- 
tations with the quantum numbers of p- meson created by the vector current. In fact, we are 
interested only in the asymptotic behaviour of these widths and masses. We shall calculate 
masses of these excitations using the Veneziano model. The relevant radial excitations all 
have the same spin and lie on the daugter trajectories to the p meson trajectory (see Fig. 
4): 

Here a' ~ 1 GeV~ 2 . 

Although there exists a working model for the resonance masses, that is in agreement 
with the experimentally observed Regge behaviour, virtually nothing is known about the 
dependence of the resonance widths on the number n of the excitation. Based on our 
experience for 2d QCD it is natural to assume the power like behaviour: 

r„ ~ ^(-n) 1/2+7 - (4.7) 

For N c — > 00 the widths r n — > 0. The constant A is some unknown proportionality constant. 
The parameter 7 parametrises our ignorance of the 4d resonance widths. The only known 
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theoretical estimate of 7 was made by Green and Veneziano [^G] on the basis of the string 
model, and their result is 7 = —1/2. Here we shall parametrise the widths by the apriory 
arbitrary parameter 7. Using eq. f£7| one immediately obtains the q 2 dependence of the self 
energy part of the inverse meson propagator: 

ImE(g 2 + ze) = A g 2(i+7). ( 4 _g) 

Then one can use the dispersion relations to obtain S(g 2 + ie). In order to do it we shall 
need to take the integral 



N c tc J so s' - q 2 

. As usual we calculate this integral for negative q 2 and then analytically continue to positive 
q 2 . Making substitution s'/(—q 2 ) —v/(l — v) we obtain that 

I = A{ ~f T ){1+7) (B(j, 1 - 7) - -(*/(! + x)yF(j, 7) 1 + 7, x/(l + *)). (4.9) 

J\ C TT 7 

Here x = —so/q 2 , and F(a,b,c, z) is the hypergeometric function. Then for large q 2 and 
noninteger 7 one obtains: 

S(g 2 ) ~ - A (~g) 2(1+7) . 71 = g 2 (^)(ctg(7r 7 ) + i)(A/N c ). (4.10) 
I\ C 7T sm(7T7j 

Note that for the 7^0 one obtains 

S(g 2 ) = ^-g 2 log(-g 2 Ao). (4.11) 
Finally for integer negative 7 one needs the dispersion relation without substractions: 

^'w/t^^ 1 ^ (412) 

Note that for each case S(g 2 ) is real for q 2 < and acquires the imaginary part for q 2 > 0. 
Now one can immediately obtain the polarisation operator: 

f 2 

U{q 2 ) = v ' ■'" 



n q 2 — n/a' + S(g 2 + ie) 
= / 2 «V(-(g 2 + S(g 2 + ^))«') 

(4.13) 

Here we used the fact that, analogously to QCD2 in order for the sum over resonances to 
be consistent with the leading asymptotic behaviour, f 2 must be independent of n for large 
n. 

For q 2 < (Euclidean domain) the polarisation operator is well approximated by OPE: 
H(g 2 ) ~ /Vlog(-gV) + +0(l/g 2 , l/i\Q... (4.14) 
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The fact that we evidently get here the terms of the order 1/q 2 reminds us that some of our 
results are model dependent. It is evident, that by changing n dependence, i.e. taking into 
account that f n are constant only for large n, one can also change the behaviour of subleading 
terms in OPE. We do not write these terms explicitly due to their model dependence. 

Consider now the q 2 > case -Minkowski domain. We can analytically continue eq. |4.14 
into that domain. The polarisation operator evidently acquires the imaginary part: 

Imn(g 2 ) ~ (/V)tt + ... (4.15) 

All the functions that we obtain by analytical continuation of the OPE are evidently 
smooth. In addition, in the Minkowski kynematics one gets the oscillating part. For suffi- 
ciently large q 2 and nonpositive 7. 



Here 

for 7 = and 



Imn(g 2 ) ~ 7ra7 2 imctg(7r(gV + iy)) (4.16) 
y = Aq 2 a'/N c 



y = A\q 2 \ 1+ y(N c sm(n 1 )) (4.17) 
for noninteger 7. For integer negative 7 

V = (A ' N ^X^ (418) 

Note that for negative 7 < — 1 one gets nondecreasing oscillations. This shows that in 
real QCD one must have the constraint 

7 > -1. (4.19) 

Moreover, we also see that 

7 < 

. Indeed, we have seen that if 7 is positive, than there appears the term in E that leads to 
the rise of the real part of the polarisation operator quicker than q 2 (although this term is 
suppressed by 1/N C ). This will lead to the curving upwards of Regge trajectories that looks 
quite unlikely at present. 

Numerically, the behaviour of the polarisation operator and its imaginary parts are 
depicted in Figs. 3 for the sample value A/N c =0.3 and 7 = —1/2 (7 = behaviour is very 
similar to the one given by Fig. 2). In these figures we depict as FO the imaginary part of the 
i/j function, as Fl the ctg approximation, and as F2 the exponentially decreasing asymptotic 
oscillations. Already starting from q 2 > 1/a' the full polarisation operator is relatively well 
approximated by ctg. Starting from the scale 

q 2 > (AT c /(A«')) 1/(1+7) (4.20) 

the full polarisation operator operator is represented by the OPE plus exponentially decreas- 
ing asymptotics. For scales much larger than the one given by eq. |4.20| the local duality is 
restored and deviations from it are exponential. 
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V. CONCLUSION 



We now can answer the basic questions asked in the introduction. First, how the local 
quark- hadron duality works, and, second when one can use the OPE. 

For the 2d QCD one can treat these two questions in a model independent way (the 
deviations from the Breit-Wigner approximation can be estimated and are actually irrelevant 
for the problem at hand). We found the polarisation operator of two scalar currents and 
studied its duality properties. 

We see that in the Minkowski domain there are two characteristic scales in the problem. 
The first scale is Ai = fin, the second is given by A2 = N c nfi 2 /(2B) For the momentum 
q 2 between these two scales the oscillations dominate, and only global duality works. For 
the q 2 , larger than A2 the local duality works. The deviations from local duality are only 
exponential. Consequently, for scales q 2 >> A 2 , one can use the local duality to calculate 
physical quantaties due to the exponential decrease of the amplitude of the oscillations. 

We thus can use the OPE to calculate global averages, i.e. to calculate some integrals 
over the interval of q 2 much bigger than the distance between the resonances for almost all 
values of q 2 . However the use of OPE based on local duality is much more restricted, and 
we can say now that there is a dynamically generated scale in the problem starting from 
which the local duality is valid. 

Moreover, we have seen that in the Minkowski domain there are always corrections to 
OPE, in the best case (q 2 > A 2 these corrections decrease exponentially.) 

The next question is what can be learned for 4d QCD. Here we have seen that if Veneziano 
model and power like behaviour of the resonance widths are indeed true, the local duality 
is established for large N c in the same way as in 2d QCD. Once again, using this time 
Veneziano model, we see there two relevant scales are I /a' and (N c / (Aa')) 1 ^ 1+ ' y ' given 
by eq. |4.20| . Between these scales the spectral density is dominated by oscillations with 
decreasing amplitudes. For the q 2 larger than these scales the spectral density behaves 
like the OPE (analytically continued, term by term, as discussed in introduction) plus the 
exponentially decreasing oscillating term, that decreases like exp(— 2Aa'q 2 ( 1+J ^) where A is 
some proportionality constant. 

We can compare ourt results with the predictions of instanton vacuum model We 
see, similar to that model that there are three characteristic areas of momenta squared, 
divided by two scales. The deviations from duality (i.e. corrections to the OPE) decrease 
exponentially for the scales larger than the biggest of these two scales, like in Ref. H. 
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FIGURES 




FIG. 1. Comparision of the decay width with the trial function. The oscillating curves are 
widths T a of resonances a's in (a), and r a averaged over 20 resonances in (b), while the smooth 
curve is T a = 0.442a ' 5 . (All widths are given in the unites h/(tt 3 N c ) 
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F 




(d) (e) 
FIG. 2. Polarisation 

functions in 2d QCD. Here: (a) F = Im^VH? 2 + j^Log(q 2 /q 2 ) + iBq 2 /N c )/(fi 2 TT 2 )), (b) 

Fi Im^ctgvr^ 2 + ^logq 2 /f + iBq 2 /N c )/(tt 2 ^))), (c) 

F 2 = -^(1 + 2exp(-2 J Bg 2 /7riV c ^ 2 )cos(2g 2 /(^ 2 ))), and (d) for F 1 /F , (e) for F 2 /F . We 

take N c = 10, B = 4.4 x 10~ 3 and fi = 0.01. We use here obvious units of dimension. 
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3. 



Polarisation 



functions 



m 



4d 



QCD 

2(l+7)> 



for 7 = -0.5 and A/N c = 0.3. (a) F = Im^(-(q 2 + A/ (N c ir)q 2{1+l) ctg^) + iA/(N c Tr)q 2y± ^ v ))), 
(b) Fl - TTlmctg(ir(q 2 + A/(N c n)q 2{1+l) ctg^) + iA/(N c ir)q 2(1+ ^))), (c) 

F2 = -vr(l + exp(-2^/(iV c 7r)g 2(1+7) )cos(27r(g 2 + ^/(iY c 7r)g 2{1+7) c^(7r 7 )))), and (d) for F^Fq, 
(e) for F2/F0. We use here obvious units of dimension. 
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